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Abstract 

Recently, several advanced multi-antenna radio 
communications technologies have emerged to meet 
the increased capacity demands in wireless multi¬ 
user networks. Despite their great potential, the 
extent of these techniques’ practical applicability 
still remains questionable, since they have to face 
either backhaul limitations or cost and hardware 
constraints. In this paper, we propose a new sys¬ 
tem solution which includes network architecture, 
antenna technology and radio transmission proto¬ 
col to reduce drastically the hardware complexity 
and cost as well as the channel state information 
/ user data feedback requirements of multi-user 
multi-antenna wireless networks. We focus on the 
forward link of an interference channel in a cloud 
radio access network setup wherein an arbitrary 
number of remote radio heads are each equipped 
with a single radio frequency module parasitic an¬ 
tenna array and wish to send data to their respec¬ 
tive single-antenna user terminals, while co-existing 
in time and frequency. Base stations select cooper¬ 
atively the optimal combination of pre-determined 
beams prior transmission. Our proposed approach 
is able to achieve the aforementioned goals, while 
offering significant downlink sum-rate gains due to 
the available spatial degrees of freedom. 

Index terms — Remote Radio Head (RRH), 
Electronically Steerable Parasitic Array Radia¬ 
tor (ESPAR) Antenna, Channel State Information 
(CSI), Zero-Forcing (ZF) Precoding, Interference 
Channel (IFC). 

1 Introduction 

Over the last decade, an enormous growth of mobile 
data traffic has been witnessed. As several studies 
indicate 00 . this trend is likely to continue in the 
years to come in an even more abrupt pace. In 
view of the well-known scarcity and high cost of 
the radio spectrum, current cellular mobile broad¬ 
band networks incorporate multiple-input multiple- 
output (MIMO) technology as a response to the 


increased capacity demands. 

Traditional single-user MIMO (SU-MIMO) fo¬ 
cuses on the physical layer performance of the com¬ 
munication between a base station (BS) and a mo¬ 
bile terminal (MT) over a point-to-point link. It 
exploits the additional degrees of freedom (DoFs) 
provided by the use of multiple antennas to enhance 
spectral efficiency through spatial multiplexing of 
individual data streams or to reduce co-channel in¬ 
terference by spatially focusing transmissions and 
separating co-channel signals. 

Multi-user MIMO (MU-MIMO) is an evolution 
of SU-MIMO that enhances link-level performance 
by exploiting DoFs to spatially separate users, thus 
enabling a BS to communicate simultaneously with 
multiple mobile users (MUs) over a multipoint-to- 
point (uplink) or a point-to-multipoint (downlink) 
channel. 

Recently, new MIMO paradigms, which aim at 
improving system-level performance, have been in¬ 
troduced. Their common characteristic is that they 
increase significantly the available DoFs to boost 
the performance of the system. More specifically, 
cooperative MIMO |;3i|4] enables the cooperation 
between different BSs to mitigate inter-cell interfer¬ 
ence and increase area spectral efficiency and sys¬ 
tem capacity, while massive MIMO [5] [6 makes use 
of an excessive number of antennas at the cell site 
to orthogonalize MUs in the spatial domain. 

It is expected that both these technologies will 
be an integral part of future fifth generation (5G) 
systems. However, their transition from theory into 
implementation requires the addressing of some im¬ 
portant practical issues. 

Cooperative MIMO refers to a collection of tech¬ 
niques with varying level of cooperation, from net¬ 
work MIMO which requires the sharing of user 
data and global channel state information (CSI) 
between the individual BSs to cooperative beam¬ 
forming which requires only the exchange of pre¬ 
coding matrices between the cooperating nodes. In 
practice, backhaul latency limits the effectiveness 
of these techniques, while backhaul capacity con¬ 
straints reduce the achievable level of cooperation 
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and the corresponding performance gains. 

Massive MIMO, on the other hand, places a 
heavy burden, in terms of cost, to the mobile net¬ 
work operators due to the large number of radio¬ 
frequency (RF) chains (one per active antenna el¬ 
ement) that is required for its implementation. 
Moreover, the tight spacing of antenna elements 
may lead to reduced efficiency due to spatial corre¬ 
lations. 

In addition, as the number of BSs or/and anten¬ 
nas increases, channel estimation becomes a chal¬ 
lenging task. 

In this paper, we describe a system setup that is 
promising in its ability to address these challenges. 
We focus on the study of a radio transmission pro¬ 
tocol that takes advantage of the attributes of the 
employed technologies to achieve substantial per¬ 
formance gains, while reducing significantly system 
complexity. 

The structure of the remainder of the paper 
is as follows: In Section [2j the considered sys¬ 
tem architecture and antenna technology is pre¬ 
sented. Next, the system and channel model 
is described in Section [3] We continue in Sec¬ 
tion [3] with a presentation of the proposed ra¬ 
dio protocol. In Section [5] we present numeri¬ 
cal performance results for different flavors of the 
proposed transmission scheme and compare them 
against performance bounds. We then discuss in 
Section [6] about the complexity reduction that is 
achieved with the considered radio protocol as well 
as about the relevant complexity-feedback over¬ 
head and complexity-performance tradeoffs. In 
Section [7] we study the performance of the trans¬ 
mission scheme when CSI is imperfect. Finally, in 
Section [8] we present our conclusions and discuss 
about future extensions of our work. 


2 System Architecture and 
Technologies 

According to the discussion in Section [lj we define 
three main criteria that are important to be met in 
order to enable the use of Cooperative / Massive 
MIMO in practice: 

1. Low-cost implementation of systems with 
many antennas. 

2. Low-delay communication between cooperat¬ 
ing nodes. 

3. Minimization of the information that has to be 
exchanged between the cooperating base sta¬ 
tions. 

In order to fulfill these requirements, all aspects 
of system design have to be taken into account, that 
is, network architecture, antenna technologies, and 


radio transmission protocols. In this Section, we 
describe two proposed components that can meet 
the first two of these objectives. Next, we will focus 
on a proposed radio scheme that can be employed 
by this system to reduce feedback complexity. 


2.1 ESPAR Antennas 

Advanced parasitic antenna arrays, often called 
electronically steerable parasitic array radiators 
(ESPARs), provide multi-antenna functionality us¬ 
ing fewer active elements (and, as a consequence, 
fewer RF units) than conventional antenna arrays. 
They accomplish this by adjusting the loads of par¬ 
asitic elements that are placed in the vicinity of the 
active element (s), thus controlling the correspond¬ 
ing induced currents caused by mutual coupling to 
form and steer beams |7 . This concept is illus¬ 
trated in Fig. [T] 

ESPAR technology can be viewed as an at¬ 
tractive enabler of massive MIMO, since it allows 
the implementation of large antenna systems with 
much reduced complexity, size and cost. 


2.2 Cloud-RAN/RRH 

Cloud Radio Access Network (Cloud-RAN or C- 
RAN) [s] |91 is a split base station architecture 
which separates the baseband unit (BBU) from the 
radio unit. More specifically, in C-RAN, remote 
radio units (RRUs), also called remote radio heads 
(RRHs), remain at the cell site, while BBUs are 
centralized and virtualized using cloud technology. 
Baseband and radio resources are allocated to “vir¬ 
tual base stations” in real time according to pro¬ 
cessing and radio coverage needs. Remote radio 
heads are connected with the centralized BBU pool 
through optical fibers. Protocols such as Com¬ 
mon Public Radio Interface (CPRI) |[10] are used 
to enable the communication between the virtual¬ 
ized baseband pool and the RRHs over this new 
network segment, which is often referred to as the 
mobile fronthaul. 

C-RAN is hence considered to be a com¬ 
pelling candidate technology for bringing coopera¬ 
tive MIMO communication into reality, due to the 
centralization and virtualization of the baseband 
and the use of optical transmission technologies 
at the fronthaul which have the potential to meet 
the stringent delay and bandwidth requirements of 


these advanced techniques 11 . 


3 System and Channel Model 

In this Section, we introduce the system and chan¬ 
nel model of our setup. Before getting in the de¬ 
tails, let us present the notation that we follow 
throughout the paper. 
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Figure 1: 


The ESPAR paradigm for beam shaping in the analog domain with digital control 
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Notation: With a, a, and A we denote a (in 
general, complex-valued) scalar, vector, and ma¬ 
trix, respectively. Arepresents the (i,j) (row, 
column) element of A. The Hermitian (transpose 
conjugate), determinant, and rank of A are de¬ 
noted by At, det(A) and rank(A), respectively, 
while the inverse and the Moore-Penrose pseudo¬ 
inverse of A are denoted by A' 1 and A + , respec¬ 
tively. Tr(A) represents the trace of A and ||A||j? 
denotes its Frobenius norm, whereas diag(a) repre¬ 
sents a diagonal matrix with the elements of vec¬ 
tor a at its main diagonal. |ja|j stands for the Eu¬ 
clidean norm of vector a. 0 ra and I n denote the 
n x n zero and identity matrix, respectively (the 
subscript indicating the dimension of the matrix 
may be omitted whenever is irrelevant). C de¬ 
notes the set of complex numbers, E{-} represents 
the expectation operator, X ~ CM (0, a 2 ) refers to 
a complex-valued random variable (RV) X follow¬ 
ing the Gaussian distribution with zero mean (i.e., 
E{X} = 0) and variance a 2 , and a ~ CM (0, a 2 T) 
represents a complex-valued Gaussian vector with 
mean matrix E{a} = 0 and covariance matrix 
R aa = E {aal} = a 2 1, i.e., the elements of a are 
complex-valued RVs with zero mean and variance 


3.1 System Model 

We consider the downlink of a C-RAN system with 
K RRHs, each representing an individual BS and 
being equipped with a single-RF ESPAR antenna, 
where each one of them wishes to communicate 
with a user equipment (UE) having a single an¬ 
tenna. All transmissions take place simultaneously 
and over the same frequency band. 

Note that due to the C-RAN architecture, the 
BSs can efficiently share information. The MUs, 
on the other hand, do not cooperate. 

Each ESPAR antenna is able to generate L dis¬ 
tinct predetermined beams. According to some cri¬ 
terion, the best A'-tuple of beams is selected for 
transmission from the L K different beam combina¬ 


tions in total, with each one of these beams being 
generated at a different RRH, as seen in Fig. [2] 
(The beam selection criteria are presented in Sec¬ 
tion |4j) 

The corresponding “beam-domain transformed” 
channel between these K transmitter-receiver (TX- 
RX) pairs is equivalent to a single-input single¬ 
output (SISO) A'-user interference channel (IFC), 
since each TX has a single active antenna. How¬ 
ever, we should not ignore the fact that the TXs are 
equipped with parasitic antenna arrays and make 
use of beamforming, in contrast to conventional 
SISO TXs which utilize omni-directional antennas. 

Mathematically, this SISO IFC with K TXs and 
K RXs having {N T , k } k=l = 1 and {Nn.k} k= : = 1 
active antennas each, respectively, can be viewed 
as a multi-user multiple-input single-output (MU- 
MISO) A-user IFC formed by a composite TX with 
Nt = l Nr,k = K antennas and K RXs, each 
with {Nft. t h} k=1 = 1 antenna. 

We assume narrowband transmission (i.e., flat¬ 
fading channel), such that no inter-symbol inter¬ 
ference is caused by multipath. This is the case, 
for example, in indoor environments with relatively 
small delay spread as well as when orthogonal fre¬ 
quency division multiplexing (OFDM) transmission 
is employed to convert a frequency-selective chan¬ 
nel into parallel frequency-flat sub-channels. We 
also assume block-fading, such that the channel re¬ 
mains fixed during the transmission of each symbol. 
Therefore, we can omit time dependence in the sys¬ 
tem model. 

Under these assumptions, the baseband received 
signal at user k, where k = 1,..., K, for a selected 
beam combination , can be expressed as 

K 

yk = ^2 h k , m Sm+n k , k = 1,... ,K (1) 
m= 1 

where s m £ C is the complex symbol transmitted 
by Tx m , y k £ C is the complex received symbol at 
Rxfc, n k ~ CM(0, cr 2 ) is complex additive Gaussian 
noise with zero mean and variance af, and h k , m is 
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Figure 2: System model. Note the transformation of the channel to a “beam-domain” interference / 
broadcast channel. 


the complex channel coefficient from Tx m to Rx*,. 
Equivalently, we can write 

y k = hls + n k , k = l,...,K (2) 

where hfc £ C Kxl is the vector of the channel coef¬ 
ficients between Rx*, and {Tx m }^ =1 and s £ C Kxl 
is the vector of transmitted symbols. 

The input-output model of the overall channel is 

y = Hs + n, (3) 

where y £ C i ^ xl is the vector of received sym¬ 
bols; n £ C Kxl is a complex-valued Gaussian 
noise vector with n ~ CAf(0, cj 2 Ik), i.e., with 
mean matrix E{n} = 0/,- and covariance matrix 
R„„ = E{nnf} = (t 2 Ik; and H £ C KxK is the 
channel matrix expressed as 

H=[h 1 , ..., h K ] f . (4) 

The total transmitted power is constrained to P, 
i.e., 

Tr(R ss ) =Tr(E{sst}) < P, (5) 

where R ss = Ejssl} is the covariance matrix of 
the transmitted signal vector s. 

For convenience, and without loss of generality, 
we assume that all transmitted symbols {sfcjfcLi 
have equal variance a 2 s = P/K = 1 such that R ss = 
of I k = I k and P = Tr (R ss ) = K. 

3.2 Channel Model 

In order to capture the effect of the beam radiation 
patterns, we consider, without loss of generality, 
a geometry-based single-bounce scattering statisti¬ 
cal model [l2] [l3]. In the single-bounce approach, 
each transmit/receive path is broken into two sub¬ 
paths: transmitter-to-scatterer and scatterer-to- 
receiver (described by their direction of departure, 
direction of arrival, and path distance). The scat- 
terer itself is modeled typically via the introduction 


of a random attenuation and phase shift. In this 
paper, we assume that the scatterers are randomly 
distributed on the surface of a sphere that is placed 
in the middle of the distance between the TXs and 
the RXs, as seen in Fig. [3] This model has been se¬ 
lected because, while it is relatively simple in com¬ 
parison with other models, it captures small-scale 
fading and shadowing effects sufficiently well. 

In order to facilitate performance analysis, we 
normalize the channel matrix. More specifically, for 
fixed H (e.g., when a given channel realization over 
the coherence time of the channel is considered), 

the normalized channel matrix H is given by 14 


H = aH, 


( 6 ) 


where, in the approach that we followed, the nor¬ 
malization constant is expressed as 


2 = N t N r = K 2 
" l|H||* ||H|| 2 f ’ 


(7) 


such that the power gain of the channel is 


Tr 




NtN r = A 2 . 


( 8 ) 


Of course, when the channel coefficients are RVs 
(e.g., longer time intervals are of interest), the ex¬ 
pectation over H should be taken in the appropriate 
expressions. 

Note that in Eqs. 0 0 ’ as well as throughout 
the paper, we denote the normalized channel matrix 
simply as H, for convenience. 


4 Radio Transmission Proto¬ 
col 

System operation is divided into three phases, 
namely, the learning phase, the selection phase, and 
the transmission phase. In the learning phase, each 
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Figure 3: Single-bounce scattering model with 100 scatterers for a system with K = 2 and L = 4 beams 
per transmitter. 


UE acquires information about the quality or the 
coefficients of each channel formed by the L K dif¬ 
ferent beam combinations at the TXs. In the se¬ 
lection phase, the BSs select jointly the best K- 
tuple of beams, according to the information that 
has been fed back to them by the UEs. Finally, in 
the transmission phase, communication takes place. 
The transmission strategy depends on the channel 
information available at the TXs. 

4.1 Learning Phase 

There exist two different types of channel informa¬ 
tion that a UE may send back to the corresponding 
BS: 

1. Its signal to interference plus noise ratio 
(SINR) for each beam combination at the TXs. 

2. The CSI for the resulting “beam-domain trans¬ 
formed channel”, i.e., the complex-valued co¬ 
efficients for the self- and cross-channels of the 
total K x K channel. 

4.2 Beam Selection Phase 

According to the type of channel information that 
the BSs acquired through feedback, the following 
beam selection criteria apply: 


We should note that, even though our goal in 
this work is to generate beams using parasitic an¬ 
tenna arrays, since ESPAR technology reduces the 
cost, size, and complexity of the system, the general 
principle holds also for beams generated by conven¬ 
tional antenna arrays. 

4.3 Transmission Phase 

In this phase, the transmission strategy that will 
be followed to send information symbols over the 
selected beams is determined. The performance 
metric that we wish to increase as much as possi¬ 
ble while maintaining complexity and CSI feedback 
overhead low is the sum-rate throughput. 

If beam selection is based on the measured SINRs 
at the RXs (i.e., the beams that maximize the sum 
of the received SINRs have been selected), then 
there is no additional processing taking place at 
the TXs. The sum-rate in this non-precoded trans¬ 
mission case is given by 


K K 


Rnp 

=s s r ) =x> 2 (i+7r ) ). 

(9) 


k =1 k =1 


where 

flr=iog 2 (i+7r ) ) 

( 10 ) 

is the rate achieved by user k and 7 k is the SINR 


of user k which is expressed as 


Beam Selection Rule #1: Select the K- 
tuple of beams that results in the largest sum of 
SINRs, i.e., the beam combination that maximizes 
2^k= i Tfc- 


Beam Selection Rule #2: Select the K- 
tuple of beams that results in the equivalent 
channel matrix with the largest product HHd. 


(NP) 

Tk 


\ h k,k\ 2 

Tnn^k \ h k,m \ 2 ^ 

_ M _ 

Srra^fc I hk,m | + &n 


( 11 ) 


since {ct 2 }^ = i = CT 2 = l. 

If, on the other hand, the channel matrices have 
been fed back to the TXs (i.e., the beams that result 
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in an equivalent channel with the largest product 
HHt have been selected), then the TXs jointly pre¬ 
code the transmit vector s, thus transforming the 
MU-MISO IFC into a MU-MISO broadcast channel 
(BC). 

We assume that linear precoding is employed by 
the system, in order to avoid the high RX complex¬ 
ity required by more advanced transmission tech- 
Linear precoding includes a family of 
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mques 

simple but sub-optinral pre-processing techniques 
that exploit CSI at the transmitter (CSIT) to im¬ 
prove MU-MIMO performance-e.g., to increase the 
sum-rate throughput or to minimize the aggregated 
bit error rate. (Unless stated otherwise, perfect 
CSIT is assumed.) 

A linear precoder generates a precoded signal 
vector as a linear transformation of the original 
symbol vector: 

s' = Ws. (12) 

Hence, the received signal vector is expressed as 


y = Hs' + n. 

That is, 

y = H W s + n . 

[Kx 1] [K x K] [KxK] [K x 1] [ifxl] 

Equivalently, we can write 

K 

y k = hj,w fe Sfe + ^2 h[.w m s m + n k , 


(13) 

(14) 


(15) 


where h*, £ C Kxl , w k £ C Kxl and s k are the 
channel vector, precoding vector, and data stream 
of user fc, respectively; H is the channel matrix; and 
W is the precoding matrix expressed as 

W = [wi • • • w k\ ■ (16) 


Note that the second term at the right-hand side 
of Eq. (15) represents the multi-user interference 
(MUI). 

The SINR of user k for this linear precoding 
scheme is given by 


7i LP) 


E 




m^k 


hfcWji 


(17) 


The expressions of the fcth user’s rate and 

the sum-rate of the system J?lp are similar with 
the ones given in Eq. (101 and ([9]), respectively, for 
the non-precoded system. 

The Tx power constraint is 


Tr(E |s' (s') f }) =Tr(fi{(Ws) (Ws) t j) 

= Tr(E{Wss t W t }) 

= Tr(WE{ss t }W t ) 

= Tr (WR ss W t ) < P. (18) 


We propose the use of zero-forcing (ZF) precod¬ 
ing, which provides a promising compromise be¬ 
tween complexity and performance 15), especially 


at the high signal to noise ratio (SNR) regime 16 


The reasoning behind ZF precoding is to employ 
a linear transformation to the transmit signal vec¬ 


tor, according to Eq. (12), such that MUI becomes 


null, i.e., so that each user k = 1,..., K receives no 
interference from the signals intended for the other 
users 
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hlw< ZF ) 

k m 


= 0, m/l 


(19) 


Interference suppression at the TX is important in 
this setup, since single-antenna, non-cooperating 
RXs are unable to eliminate interference-the best 
they can do is to treat it as noise. 

A common approach to accomplish this goal in¬ 
volves the inversion of the channel matrix at the TX 
in order to create orthogonal channels between the 
TX and the RXs (i.e., to diagonalize the effective 
channel). Channel inversion implies, under the flat¬ 
fading assumption, the multiplication of the origi¬ 
nal transmit vector signal with the right Moore- 
Penrose pseudo-inverse of the channel matrix H, 
such that the matrix of the composite channel is 
the identity matrix. That is, 


W ZF = H+ => 


s' = W ZF S = H+s 

HW zf = HH+ = I 


( 20 ) 


where 

H+ = H t (HH t ) _1 . (21) 


In practice, though, we typically normalize the 
precoding matrix W zf in order to set the transmit 
power (after precoding) to a fixed value, indepen¬ 
dent of the channel H, according to the given power 
constraint. An often used normalization, which is 
referred to as equal receive power (ERP) normal¬ 
ization, takes the following form 17 : 


W ER p = V^H+ = (HHt) \ (22) 


Under our assumption a'l = 1 => E {sst} = R ss = 
I k => Tr (R ss ) = K and by forcing the total 
transmit power after precoding to remain equal to 
P = K, that is, 


Tr (WerpR^W^p) = K, (23) 

the power normalization factor (3 is given by 118; 


Tr ((HHt) 

-) 

1 


(HHt) -1 

2 

1 



Ef=iV(A l) 


(24a) 

(24b) 

(24c) 
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where A& is the fcth singular value of H. 

When the channel matrix is square, as it is in our 
case where N T = K , then the precoding matrix is 
simply given by 19]: 

W ERP = (25) 


with the norm of the corresponding column vec¬ 
tor [20 


Wrtp — 


F(:> k) 


|F(:,fc)||’ 
The SINR in this case is 


k = 1,..., K. 


(31) 


In any case, the corresponding received signal 
vector is given by: 


(ETP) 

Tk 


1 

a n IIF(:, fc)|| 2 


(32) 


y = Hs' + n 
= H (W ERP s) + n 
= (HW E rp) s + n 
= \J~j3s + n. (26) 


Note that the scaling factor in Eq. (26) may lead 
in reduced SNR at the receiver if the channel matrix 
is ill-conditioned, i.e., if one of its singular values is 
very large. 

The SINR at user k is given by Eq. ( |l7| ) by set¬ 
ting the MUI at the denominator equal to zero and 
noting that HWerp = y/(3T- 


% 


(ERP) _ P 


(27) 


Therefore, all MUs experience the same SNR as 
expected, since the scaling factor is the same for all 
transmitted signals, and achieve the same rate 

4 ERPI =log 2 (l+7f HP >)' < 2s > 


and the sum-rate is given as 

-Retp = 53 log 2 (l + 7fe ETP ' ) 
fc=1 

K / 


= 53 lo s2 1 + 


*;=1 


^l|F(:,*)ir 


(33) 


Both these normalization methods set MUI to zero 
and have the same total transmit power. These ZF 
precoding design schemes aim at minimizing trans¬ 
mit power. Another design approach is to set 

W WF = Ht (HHt)^ 1 diag (y^, • • •, v 7 ^) 

(34) 

and perform water-filling (WF) to optimally allo¬ 
cate transmit power per antenna {P& }* =1 such that 
the sum-rate is maximized, under the total power 
constraint 

K 

53l|w fe ||P fc = P. (35) 

k =1 


Hence, the sum-rate of this ZF precoding scheme 
is given by: 


Repp 


e*t p) 

k =1 

£M 1 + ^ erp ’) 

5>4 + £) 

fc=1 ' "' 


K log 2 
K log 2 



P 



1 

P 

1 



(29a) 

(29b) 


Another commonly used normalization, which re¬ 
sults in higher sum-rate than ERP normalization, 
is the so-called equal transmit power (ETP) nor¬ 
malization which is obtained by setting 

F = H + =H t (HH f )" 1 (30) 

(or F = H 1 in case of a square channel matrix) 
and then dividing the elements of column k of F 


However, this optimal ZF method is not applica¬ 
ble in the considered system setup due to the fact 
that it does not only increase the complexity of the 
cooperation between the BSs, but it requires also 
cooperation between the UEs (i.e., joint decoding). 


5 Performance Results 

In this section we evaluate the sum-rate perfor¬ 
mance of the proposed radio protocol via numer¬ 
ical simulations for the system setup described in 
Section [3] with K = 2 and L = 4. 

The performance results in Fig. [4] represent er- 
godic sum-rates achieved over a range of target 
SNR values at the receiver, from OdB up to 30dB. 
These results have been obtained after 1,000 simu¬ 
lation runs by taking the expectation of the corre¬ 
sponding sum-rate equations. Also, in each simula¬ 
tion run, 100 sub-runs are used for the normaliza¬ 
tion of the channel matrices, as described in Sec¬ 
tion [372] 

We consider transmission over the “beam- 
domain” channel with and without ZF precoding. 
That is, we study non-precoded transmission over 
the beams (which have been selected according to 
the SINRs that the MTs have fed back to the BSs) 
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Figure 4: Sum-rate of different transmission schemes for the case of K = 2 interfering links. 


as well as the case where the MTs feed back the 
complex coefficients of the equivalent 2x2 “beam- 
domain” channel and ZF precoding is incorporated 
in the transmission process. 

For comparison purposes, we illustrate also the 
performance of the system when communication 
does not take place over a selected pair of beams, 
but we have instead conventional (non-precoded or 
precoded) “omni-channel” transmission. Of course, 
this type of non-precoded transmission would not 
take place in practice, since it is not a multi-user 
communication scheme, i.e., it does not take into 
account MUI; it is just included in the simulation 
as a lower performance bound. 

In our simulations we have set, for mathematical 
tractability and without loss of generality, P = K. 

We note that non-precoded transmission over the 
“beam-channels” outperforms significantly “omni¬ 
channel” transmission methods, even when the lat¬ 
ter employ ZF precoding, especially at high SNR 
values. This is due to: (a) The gain that is intro¬ 
duced to the system as a result of the use of beams. 
(b) The fact that the TXs are informed about the 
SINR at the RXs and select the best beam combi¬ 
nation each time. 

We also note that the knowledge of the CSI at 
the transmitter for the equivalent “beam-domain 
transformed” channel and the joint precoding of 
the transmit vector that is performed based on that 
CSIT further improves the performance of the sys¬ 
tem at the high SNR regime. 

Finally, it is worth noting that the sum-rate of 
the transmission schemes that incorporate ZF pre¬ 
coding keeps increasing linearly within the consid¬ 
ered SNR range, whereas the non-precoded tech¬ 
niques experience an expected flooring of their per¬ 
formance due to the residual interference that they 
incur, which is not accounted for. 


6 Complexity-Overhead and 
Complexity-Performance 
Tradeoffs 

In Section [2j a brief overview of ESPAR and C- 
RAN technologies was given. In this Section, we 
present the complexity reduction accomplished by 
the use of the proposed radio transmission scheme 
described in Section [4j mainly in terms of channel 
estimation / channel information feedback, and the 
relevant complexity-overhead as well as complexity- 
performance tradeoffs. But first, let us summarize 
the main points of Section [2] 

• ESPAR antennas reduce the size, complexity 
and cost of antenna arrays due to the fact that 
they use fewer active antenna elements and RF 
units than conventional antenna solutions and 
thus can be viewed as an enabler of massive 
MIMO. 

• C-RAN facilitates the efficient cooperation be¬ 
tween individual BSs in Coordinated Multi- 
Point (CoMP) setups since the centralization 
and virtualization of the BBUs as well as the 
use of optical fibers at the fronthaul for the 
connection of the BBUs with the RRHs is ca¬ 
pable of meeting the stringent delay and band¬ 
width requirements that are imposed by this 
family of cooperative MIMO communication 
techniques. 

The interested reader is encouraged to refer to 
the corresponding references cited in Section [2] for 
a more detailed presentation of the benefits and 
drawbacks / challenges of these technologies. 

As we mentioned in Section [3j our focus is on 
the downlink transmission in a system setup with 
K cloud-based RRHs/BSs, each equipped with 
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a single-RF ESPAR, and K single-antenna UEs. 
Since the K RRHs can be viewed as a composite 
BS with K antennas due to the C-RAN architec¬ 
ture, MU-MISO communication techniques over a 
BC are applicable. 

The use of ESPARs instead of conventional an¬ 
tenna arrays reduces significantly channel estima¬ 
tion complexity due to the decreased number of ac¬ 
tive antenna elements which is translated into a cor¬ 
respondingly decreased number of direct- and cross¬ 
channels. In the remaining Section, we will com¬ 
pare “omni-channel” and “beam-channel” trans¬ 
mission methods in terms of complexity assum¬ 
ing in both cases the use of ESPARs at the BSs, 
but we should keep in mind this important note, 
that is, that the replacement of antenna arrays by 
ESPARs already provides complexity reduction re¬ 
garding channel estimation. 

Initially, let us assume conventional transmission 
over “omni-channels”. The optimum precoding 
technique in terms of sum-rate is dirty paper cod¬ 
ing (DPC) [l3j. However, the computational com¬ 
plexity of this non-linear precoding scheme makes 
it impractical. ZF precoding, on the other hand, is 
a low-overhead scheme that completely eliminates 
MUI, but it amplifies noise power. Minimum mean 
square error (MMSE) considers noise in the precod¬ 
ing design process, leading to better performance 
than ZF in low SNR conditions, at the cost of in¬ 
creased complexity 13 . Hence, ZF precoding is the 


MU-MISO transmission method with the smaller 
overhead. 


ZF precoding requires full CSIT. Since we have K 
active antennas at the composite TX (and recalling 
that we are considering narrowband channels), the 
equivalent channel is described by a K x K matrix. 
In other words, each one of the K UEs should feed 
back to its corresponding BS a vector of size K 
containing the coefficient of the direct channel with 
that BS and the K — 1 coefficients of the cross¬ 
channels with the other AT — 1 BSs. 


If we assume non-precoded transmission, as a low 
complexity bound, then each UE should feed back 
to its BS only its SINR. Thus, the composite BS 
would have to collect AT SINR values (real num¬ 
bers), as opposed to the K 1 2 channel coefficients 
(complex numbers) required in the ZF precoding 
scenario. Of course, as we mentioned in the previ¬ 
ous Section, such a transmission scheme would not 
be used in practice due to its poor performance. 

Now, let us consider the proposed radio protocol. 
Each one of the AT RRHs is able to generate L dis¬ 
tinct beams. Thus, there exist L K possible beam 
combinations. Assuming non-precoded transmis¬ 
sion (i.e., the beams to be used during transmission 
are selected according to criterion #1 described in 
Section [4]), each UE should feed back to its BS its 
SINR for each one of the beam combinations. That 


is, each one of the K UEs should send to its corre¬ 
sponding BS L K SINR values. Therefore, K ( L K ) 
real values will be sent back to the composite BS in 
total. In contrast, when ZF precoding is used over 
“omni-channels”, K 2 complex values have to be 
send back to the composite BS. However, in the lat¬ 
ter case, accurate channel estimation is challenging 
when the number of BSs (and active TX antenna 
elements) is high, while the measurement of the re¬ 
ceived SINR value at each UE is a rather trivial 
task. Also, while non-precoded, SINR-based trans¬ 
mission over “omni-channels” presents the lowest 
overhead (the feedback of only K real values is re¬ 
quired), it is not suitable for multi-user communi¬ 
cation as we have already mentioned. 

Let us turn our attention into ZF precoding- 
based transmission over the “beam-channels” (i.e., 
beam selection criterion is applied). For each 
beam-combination, the required CSIT is repre¬ 
sented by a if x A" complex-valued matrix. Thus, 
in total A' 2 (A a ) complex values have to be fed 
back to the composite BS instead of only K 2 com¬ 
plex values that are required in the corresponding 
“omni-channels” case. However, a performance- 
complexity tradeoff is expected due to the “beam- 
gains” and the selection of the optimal “compos¬ 
ite beam-domain channel”, as we have seen in Sec¬ 
tion [5] 

Let us summarize: Non-precoded and precoded 
transmission over “beam-channels” results in an in¬ 
crease of the feedback overhead by L K in com¬ 
parison with the corresponding “omni-channel” 
schemes due to the fact that each UE should send 
back to its BS its SINR or channel matrix for each 
beam combination. However, it outperforms trans¬ 
mission over “omni-channels” and it simplifies (or 
even eliminates the need for) channel estimation. 
More specifically: 

• “Beam-channel” communication techniques 
present a significant performance gain against 
their “omni-channel” counterparts. 

• Even non-precoded “beam-channel” transmis¬ 
sion outperforms ZF “omni-channel” transmis¬ 
sion. Thus, by using predetermined beams, 
we could simply collect SINR values instead of 
estimating channel matrices, which is a much 
simpler procedure - and therefore, it is easier 
to be used in practice, even in large setups^ 

• The increase in the channel information over¬ 
head caused by the use of beams can be partly 
compensated by the use of ESPARs instead 
of conventional antenna arrays, depending on 

1 In high mobility environments, it is possible to update 

the SINR value corresponding to each beam combination 
every time it is used to serve a user. Then, a record of time- 
windowed SINR values can be kept, thus further reducing 
the feedback requirements [21](22 23]. 
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the ratio “reduction in the number of active 
antenna elements / increase in the number of 
beams”. 


7 ZF Precoding-based Trans¬ 
mission over “Beam- 
Channels” with Imperfect 
CSIT 

So far, we have assumed perfect CSIT. In practice, 
though, various sources of error (e.g., channel es¬ 
timation errors, channel quantization errors, feed¬ 
back errors etc.) may result in imperfect CSIT. 

In this case, the imperfect channel matrix that is 
fed back to the TX is expressed as 

H e = H + E, (36) 

where H is the actual channel matrix and E is an 
additive error matrix whose entries are independent 
and identically distributed (i.i.d.) RVs that are in¬ 
dependent of H and follow a [E. y ] ~ CJ\f (0, a^) 
distribution. 

Then, the precoded signal vector is expressed as 
s' = v^W zf H+, (37) 

where the power normalization factor j3 is given by 



1 


(h b h£) 

V ' F 

1 

■£f=il/(A 2 fc )’ 

and A fc is the Ath singular value of H e . 

The received signal vector is expressed as 


(38a) 

(38b) 


(38c) 


y = Hs' + n 


= V^HH^s + n 
= \fp (H e - E) H£s + n 
= y/ps- v^EH^s + n. (39) 

The SINR of each user, given that erf = 1, is 

(ZF) _ P _ _1_ 

^ ~ P °l (<?* + 1 IP) £f=l 1/Al ’ 

(40) 


while the sum-rate is expressed as 


Rzf 


K log 2 



1 

K 2 + i/P)ELi/A 2 t 



Note that the variance of the error, er 3 , limits the 
sum-rate, i.e., as P —> oo, 


Rzf = K log 2 ( 1 + * ) ■ (42) 

We conclude this Section by evaluating the sum- 
rate performance of ZF precoding transmission over 
pre-determined beams for the considered setup in 
the case where CSIT is imperfect and for a range of 
error variance values P e from 10 -3 up to 10 _1 via 
simulations and we compare these results with the 
performance obtained when CSIT was assumed to 
be perfect. 

As we see in Fig. [5j when the error variance is 
small, the degradation in the sum-rate throughput 
is negligible and only noticeable at the high SNR 
regime. However, as the error variance increases, 
the sum-rate is decreased and at some point floors 
(the higher the error variance, the sooner this floor¬ 
ing takes place). 


8 Conclusions 

In this paper, we proposed a new radio transmission 
protocol for a system setup that could be an enabler 
of next-generation MU-MIMO systems. The pro¬ 
posed scheme performs extremely well in terms of 
the achieved sum-rate while it reduces significantly 
the complexity of the system. In the future, we 
plan to extend this work by studying larger setups 
and cases with imperfect CSI feedback. 
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